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Abstract The goal of relativistic heavy-ion collisions is to determine the phase boundary of quantum 
chromodynamics (QCD) phase transitions. Critically sensitive observables are suggested to be higher-order 
cumulants of conserved charges. The non-monotonous behavior of higher cumulants was observed at the relativistic 
heavy-ion collider (RHIC). However, it remains unclear whether these non-monotonous behaviors are critically 
related. We studied the influences of non-critical fluctuations, finite system size, and limited evolution time to 
determine if they cause non-monotonous behavior. First, we examined the minimum statistics required for measuring 
the fourth cumulant. The minimum statistic obtained using the centrality bin width correction (CBWC) method was 
25 M. We suggest using a 0.1% centrality bin in the CBWC method instead of each Nch. With a 0.1 centrality bin 
width, 1 M statistics are sufficient. We then pointed out the statistical fluctuations from the limited number of final 
particles. By assuming the independent emission of each positive (or negative) charged particle, the statistical 
fluctuations of positive (or negative) charged particles were presented by a Poisson distribution, and the statistical 
fluctuations of net-charged particles were their evolution. The obtained statistical fluctuations for net protons, net 
electronic charges, and net baryons were consistent with those from the Hadron Resonance Gas model. In addition, 
the measured cumulants at RHIC/STAR are dominated by these Poisson-like statistical fluctuations. At the end of this 
section, we suggest the pooling method of mixed events and demonstrate that the sample of mixed events accurately 
presents the contributions of the background. Dynamic cumulants were defined as the cumulant of the original 
sample minus that of the mixed sample. Dynamical cumulants were shown to simultaneously reduce the influence of 
the statistical fluctuations, centrality bin width effects, and detector efficiency. Second, because the system is finite, 
the correlation length at the critical point is not developed to infinity in contrast to the system at thermal limits. Using 
a Monte Carlo simulation of the three-dimensional three-state Potts model, we demonstrated the fluctuations of the 
second- and fourth-order generalized susceptibilities near the temperatures of the external fields of the first-, second-, 
and crossover regions. Both the second- and fourth-order susceptibilities showed similar peak-like and oscillation- 
like fluctuations in the three regions. Therefore, non-monotonic fluctuations are associated with the second-order 
phase transition and the first-order phase and crossover in a finite-size system. The exponent of finite-size scaling 
(FSS) characterizes the order of transitions or crossover. To determine the parameters of the phase transition using the 
FSS, we studied the behavior of a fixed point in the FSS. To quantify the behavior of the fixed point, we define the 
width of the scaled observables of different sizes at a given temperature and scaling exponent ratio. The minimum 
width reveals the position of the fixed point in the plane of the temperature and scaling exponent ratio. The value of 
this ratio indicates the nature of the fixed point, which can be a critical, first-order phase transition line point, or 
crossover region point. To demonstrate the effectiveness of this method, we applied it to three typical samples 
produced by a three-dimensional three-state Potts model. The results show that the method is more precise and 
effective than conventional methods. Possible applications of the proposed method are also discussed. Finally, 
because of the limited evolution time, some processes in relativistic heavy-ion collisions may not reach thermal 
equilibrium. To estimate the influence of the nonequilibrium evolution, we used the three-dimensional Ising model 
with the Metropolis algorithm to study the evolution from nonequilibrium to equilibrium on the phase boundary. The 


order parameter exponentially approaches its equilibrium value, as suggested by the Langevin equation. The average 
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relaxation time is defined. The relaxation time is well represented by the average relaxation time, which diverges as 


the zth power of the system size at a critical temperature, similar to the relaxation time in dynamical equations. 


During nonequilibrium evolution, the third and fourth cumulants of the order parameter could be positive or negative 


depending on the observation time, which is consistent with the calculations of dynamical models at the crossover 


side. The nonequilibrium evolution at the crossover side lasts briefly, and its influence is weaker than that at the first- 


order phase transition line. These qualitative features are instructive for experimentally determining the critical point 


and phase boundary in quantum chromodynamics. 


Key words Relativistic heavy ion collisions, QCD phase transition, Non-critical fluctuations, Finite system size, 


Non-equilibrium 


相对 论 重 离子 碰撞 的 一 个 主要 科学 目的 是 确定 
量子 色 动 力学 (Quantum Chromo-dynamics,QCD ) 的 
相 图 。 为 实现 这 个 目标 ,人 们 设计 并 开展 了 大 量 的 
实验 。 目 前 ,在 美国 布鲁克 海 文 国家 实验 室 正在 进 
行 的 相对 论 重 离子 碰撞 能 量 扫描 计划 (Relativistic 
Heavy-Ion Collider Beam Energy Scan, RHIC-BES ) 
II ,以 及 未 来 俄罗斯 的 NICA CNuclotron-based Ion 
Collider Facility) 、 德 的 FAIR (Facility for 
Antiproton and Ion Research) 和 中 国 的 HIAF (High 
Intensity Heavy-ion Accelerator Facility) 实 验 计 划 都 
是 希望 通过 能 量 扫描 发 现 与 QCD 相 变 相关 的 观测 
信号 ,实现 在 实验 上 确定 QCD 相 变 临界 点 和 相 边 
FR), KUE, ARSE A AE QCD 相 变 相关 的 研究 将 是 
相对 论 重 离子 碰撞 的 国际 前 沿 。 

QCD 预言 在 足够 高 的 温度 和 低 重 子 化 学 势 的 
条 件 下 ,物质 会 发 生 从 强 子 气体 到 夸克 胶 子 等 离子 


效应 的 影响 。 我 们 知道 ,临界 涨 落 是 由 于 事件 内 部 
粒子 之 间 的 关联 所 致 ,其 关联 长 度 具 有 发 散 的 特征 ， 
而 非 临界 的 涨 落 不 具有 关联 长 度 发 散 的 特征 , 非 临 
界 涨 落 有 两 个 主要 来 源 :一 是 源 于 传统 的 守恒 机 制 ， 
如 共振 衰变 , 以 及 能 量 、 动 量 、 电 蓓 数 的 整体 守恒 。 
这 些 传统 效应 的 关联 长 度 收敛 且 固 定 。 相 对 于 临界 
涨 落 , 非 临界 涨 落 通常 比较 小 ,对 于 守恒 荷 高 阶 矩 的 
影响 也 是 固定 的 。 第 二 个 来 源 是 由 观察 系统 本 身 的 
局 限 性 和 特点 所 致 ,例如 :有 限 事件 数 ( 即 :统计 量 ) 
所 带 来 的 测量 不 确定 性 一 有限 粒 子 数 引 起 的 统 
计 涨 落 ”不 同 碰撞 参数 导致 初 态 系统 太 度 的 涨 
落 、 中 心 度 宽度 、 探 测 器 的 接收 效率 与 有 限 实验 探测 
相 空 间 所 给 出 的 人 为 截断 ” ”等 。 

为 扣除 传统 关联 机 制 和 系统 效应 所 带 来 的 非 临 
界 涨 落 ,人 们 尝试 提出 了 针对 各 种 效应 的 修正 方案 。 
例如 :1) 对 于 统计 涨 落 的 影响 ,通常 采用 具有 相同 粒 


一 


体 的 相 变 。 格 点 QCD 理论 计算 表明 :在 高 温 零 化 学 
势 附 近 , 从 强 子 到 夸克 胺 子 等 离子 体 的 转变 不 是 一 
个 真正 的 相 变 ,而 是 一 个 解析 的 平滑 过 渡 ”” 。 基 于 
QCD 的 唯 象 模型 预言 :在 低温 和 高 重子 化 学 势 时 ， 
强 子 到 夸克 胶 子 等 离子 体 的 相 变 是 一 级 相 变 ,一 
级 相 变 线 的 终止 点 是 QCD 的 临界 点 ,QCD 的 相 图 
由 一 级 相 变 线 和 临界 点 完全 确定 。 


子平 均 数 的 泊 松 分 布 来 模拟 "中 ;2) 在 一 个 给 定 的 
中 心 度 宽度 中 ,实际 上 混杂 了 多 个 碰撞 参数 并 不 完 
全 相同 事件 的 合 加 ,通常 高 阶 矩 随 中心 度 宽度 变化 
不 是 平滑 的 。 为 了 减少 由 此 带 来 的 影响 ,普遍 采用 
了 一 种 被 称 为 中 心 度 宽度 修正 CCentrality Bin 
Width Correction, CBWC) MW 77°"; 3) FAT we E 


根据 格 点 QCD 理论 和 有 效 势 模型 ,守恒 荷 的 高 
阶 矩 对 于 关联 长 度 特别 敏感 , 且 在 临界 点 附近 可 能 
会 出 现 非 单 调 行为 和 符号 变化 ,因此 ,QCD 相 变 临 
界 点 的 敏感 观测 量 是 末 态 守恒 荷 粒子 的 高 阶 矩 1。 
在 RHIC/BES I 上 ,人 们 已 经 观测 到 净 质 子 数 的 四 阶 


到 ,这 会 直接 影响 粒子 数 分 布 的 峰值 以 及 累积 
为 了 消除 这 种 探测 器 效率 的 影响 ,人 们 建立 了 真 
累积 矩 与 实际 测量 累积 矩 之 间 的 对 应 关系 ,利用 对 
应 关系 进行 修正 2 。 但 是 ,对 于 一 个 真实 的 样本 ,所 
有 这 些 非 临 界 效应 是 相互 联系 的 ,很 难 单独 地 将 某 


EA 
子 碰 撞 实 验 中 ,有 部 分 末 态 粒子 没有 被 探测 器 观察 
5 


和 矩 随 碰撞 能 量 的 变化 具有 非 单 调 行为 ”"” ,但 是 ， 
非 单调 行为 并 不 是 临界 涨 落 所 独 有 的 特征 。 在 一 级 
相 变 和 平滑 过 渡 时 , 相 变 相关 的 涨 落 也 会 有 非 单 调 
PENS) ;另外 ,实验 还 观测 到 在 各 种 碰撞 能 量 下 , 净 
质子 数 的 三 阶 矩 都 是 小 于 泊 松 结果 的 负 值 , 这 与 平 
衡 相 变 所 预期 的 正定 结果 不 相符 2 3。 

为 了 确定 所 观测 到 涨 落 的 临界 特征 ,大 量 理论 
和 实验 数据 分 析 工 作 致 力 于 扣除 各 种 可 能 的 非 临 界 


种 效应 消除 掉 。 一 种 更 好 的 消除 思路 是 同时 将 这 
些 非 临界 效应 进行 处 理 , 混合 事件 方法 提供 了 得 到 
这 种 背景 的 可 能 。 

本 文 将 在 $1 讨论 这 些 非 临界 效应 的 影响 和 扣 
Ro KR ,在 $1.1 讨论 有 限 事例 数 对 高 阶 窍 ,尤其 是 
四 阶 和 六 阶 矩 的 影响 “六 ;接着 在 8$1.2 介绍 对 统计 
涨 落 的 描述 ,以 及 统计 涨 落 对 高 阶 矩 实验 测量 结果 
的 影响 ;最 后 在 $1.3 给 出 有 效 混合 事件 的 池子 方法 ， 
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定义 动力 学 高 阶 矩 为 原始 样本 的 高 阶 矩 减 去 混合 样 。” 置 所 对 应 的 温度 和 标 度 指数 率 即 为 临界 温度 和 临界 
本 的 高 阶 矩 ,并 展示 动力 学 高 阶 矩 对 扣除 非 临 界 涨 ” 标 度 指数 率 。 在 重 整 化 群 理论 中 ,临界 点 为 不 稳定 
落 的 有 效 性 。 固定 点 ,一 级 相 变 线 上 存在 不 连续 的 固定 点 。 利 用 
除了 非 临 界 涨 落 以 外 ,我 们 知道 还 有 两 个 重要 相 边 界 上 的 固定 点 行为 ,可 以 精确 决定 临界 参数 。 
且 不 可 忽略 的 因素 ,首先 ,相对 论 重 离子 碰撞 所 形成 ” 本 文 将 在 $2.2 详 细 介 绍 由 相 边 界 上 观测 量 的 固定 点 
的 系统 不 是 热力 学 极限 下 的 无 穷 大 ,而 是 只 有 几 个 ” 行为 定量 确定 其 临界 参数 的 方法 ™。 

fm 的 大 小 ,因此 ,必须 考虑 有 限 系 统 尺度 对 观测 量 对 于 非 平 衡 演 化 ,目前 还 没有 成 熟 的 非 平衡 统 
的 影响 ;同时 ,系统 的 演化 时 间 也 只 有 10 fme 计 理 论 , 非 平 衡 效 应 对 相 变 以 及 临界 涨 落 的 影响 尚 
右 ,非常 有 限 汪 ,因此 系统 很 可 能 不 能 全 部 达到 热力 ”无 法 定量 估计 。 基 于 非 平衡 演化 的 动力 学 理论 也 只 
学 的 平衡 态 , 必 须 考 虑 非 平 衡 演 化 对 观测 量 的 影响 。 ”能 局 限于 对 平滑 过 渡 区 域 进行 定量 计算 ”“"”。 要 估 


这 两 个 因素 对 于 正确 决定 临界 点 、 相 边界 在 QCD 相 。” 计 非 平 衡 效应 对 临界 点 一 级 相 边界 ,以 及 相 边 界 附 

图 上 的 位 置 至 关 重 要 。 近 的 影响 ,我们 必须 寻求 新 的 描述 非 平衡 演化 的 有 
限 系统 尺度 效应 使 临界 关联 长 度 的 发 散 受 到 效 模型 和 方法 。 

限制 , 涨 落 信号 会 偏离 热力 学 极限 下 的 数值 Pr , 例 众所周知 ,平衡 相 变 是 系统 自发 对 称 性 破 缺 时 

如 ,在 热力 学 极限 下 ,磁化 率 在 一 级 相 变 处 是 一 个 5 ”所 发 生 的 物理 现象 ,具有 相同 对 称 性 的 系统 所 发 生 

函数 ,在 二 级 相 变 处 是 发 散 的 无 穷 大 。 但 是 ,对 于 有 的 相 变 属于 同一 个 普 适 类 ,同一 普 适 类 的 系统 具有 

限 尺度 系统 ,磁化 率 在 一 级 和 二 级 相 变 处 都 是 一 个 。 相同 的 临界 指数 。 在 统计 物理 中 ,Ising 模 型 是 研究 


有 限 峰值 ,并 且 其 峰值 在 相 图 上 所 对 应 的 位 置 ,如 临 。 相 变 特征 最 有 效 的 模型 之 一 , 它 很 好 地 展示 了 真实 
界 温度 也 会 偏离 其 临界 值 , 变 为 大 临界 温度 。 本 ” 相 变 系统 的 动力 学 特征 。QCD 和 三 维 Ising 模 型 的 


文 $2.1 将 利用 三 维 三 态 Potts 模型 完整 展示 有 限 尺 ”临界 点 同属 于 Z(2) 对 称 群 的 普 适 类 ,通过 选择 

度 系统 在 相 边 界 附 近 的 这 些 特征 。 适当 的 映射 关系 ,可 以 将 Ising 模 型 的 温度 和 外 场 映 
对 于 有 限 尺度 系统 ,要 精确 确定 其 临界 点 位 置 ，” 射 到 对 应 QCD 系统 的 温度 和 化 学 势 ”。 

必须 利用 观察 量 的 有 限 尺度 标 度 性 。 我 们 知道 ,对 在 临界 点 附近 , 非 平衡 演化 的 共同 特征 是 弛 玉 


besi 


于 像 磁化 率 这 样 的 观察 量 , 它 能 够 很 好 地 用 有 限 尺 寸 间 按 动力 学 指数 的 寡 次 发 散 , 动 力学 指数 依赖 于 
度 标 度 性 来 描述 。 e 具体 的 系统 参数 和 算法 。 虽 然 , 三 维 Ising 系 统 的 动 
临界 温度 的 数值 ,而 且 根 据 其 标 度 指数 的 数值 还 力学 指数 与 QCD 系统 不 完全 一 样 ,定量 的 数值 结果 
判定 相 变 的 级 数 22。 不 具有 直接 的 可 比 性 。 但 是 , 非 平衡 演化 在 整个 相 
对 于 给 定 的 系统 和 观察 量 ,其 有 限 尺 度 标 度 性 。 ”边界 的 定性 特征 ,如 : 非 平衡 效应 对 温度 (或 相 边 界 
I 具体 形式 必须 通过 观测 量 在 不 同系 统 尺度 下 随 温 。 区域 ) 的 依赖 性 、 非 平衡 效应 是 否 会 改变 序 参 量 及 其 
度 的 变化 行为 来 确定 。 为 此 ,一 般 传统 的 做 法 是 : 先 ”高 阶 累积 矩 的 符号 等 ,这 些 定性 特征 可 以 帮助 我 们 
假定 临界 特征 参数 , 即 临 界 温度 和 临界 标 度 指数 率 ， 估计 非 平衡 效应 对 QCD 系统 的 影响 。8$3.1 首先 利 
选 自 变量 为 标 度 约 化 温度 , 画 观 察 量 在 不 同系 统 尺 。 用 Metropolis 算 法 模拟 三 维 Ising 模 型 的 非 平 衡 演 
度 标 度 下 的 曲线 ,根据 曲线 在 临界 约 化 温度 处 的 重 ”化 ,分 析 在 临 界 点 附近 非 平衡 演化 的 动力 学 特征 ,给 
合 情 况 来 判断 事先 给 定 的 特征 参数 是 否 正确 中 。 显 ”出 对 非 平 衡 态 的 定量 描述 ;然后 ,$3.2 给 出 磁化 强度 
然 , 这 种 传统 方法 存在 局 限 性 ,首先 , 标 度 曲线 重合 ”各 阶 窍 随 观 察 时 间 的 变化 ;最 后 ,8$4 给 出 了 对 所 有 
情况 大 多 是 由 肉眼 来 判断 ,无 法 定量 描述 其 标 度 的 ”工作 的 总 结 。 
精确 度 。 而 且 , 在 临界 约 化 温度 附近 ,由 于 尺度 因子 -nh 
的 加 入 ,相应 的 标 度 曲 线 重合 区 域 非常 大 ,导致 对 应 工 “ 非 临 弄 涨 落 的 影响 
的 相 变 特征 参数 有 很 大 的 不 确定 范围 。 , 
如 果 将 自 变量 从 标 度 约 化 温度 换 为 温度 , 标 度 LI 统计 量 的 影响 
曲线 将 仅 在 相 变 温度 处 交 于 一 点 , 即 固定 点 。 在 偏 如 前 所 述 , 相 对 论 重 离子 碰撞 中 临界 涨 落 敏感 
ies 度 处 ,不 同 尺度 的 曲线 彼此 分 开 , 固 定点 位 A TE ae AY re RE, EN A he ON: 
=(X), C, =((dX)’), C, = (dX )*), C, = (dX )*) - Erry, 


C; = ((dX )°) - 10 (8X )*)((dX )?), Ce = (8X )*) - 10(( dX )} y + 30((dX)? y — 15((8X )*)((dX J? y C 
中 :6XY = 对-《X》。C 是 序 参 量 的 平均 值 ;C, 是 分 ” 布 的 方差 即 o?;C,、Cs 与 分 布 的 偏 度 S$S= C/C y 
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峰 度 k= Cy/C, 直接 相关 , 偏 度 和 峰 度 是 对 分 布 非 
高 斯 性 的 量化 。 

对 于 相对 论 重 离子 碰撞 中 的 大 多 数 观测 量 , 如 : 
净 质 子 数 的 平均 值 、 二 阶 矩 等 ,我 们 都 不 需要 考虑 统 
计量 的 问题 ,因为 对 于 这 些 普 通 的 观察 量 , 目前 实验 
的 统计 量 已 经 足够 保证 其 精准 度 。 但 是 ,对 于 高 阶 
JE ,尤其 是 四 阶 或 六 阶 矩 ,精准 测量 要 求 的 统计 量 将 
会 对 应 地 成 究 指 数 增长 。 因 此 ,我 们 必须 首先 回答 


计算 中 心 度 在 0%~5% 的 净 质 子 四 阶 矩 ,结果 发 现 ， 
《xo’) 随 着 统计 量 而 增加 ,违背 了 中 心 极限 定理 ,其 
原因 是 每 个 Nch 的 事件 数 太 少 , 没 有 达到 中 心 极限 
其 实 ,在 现 有 统计 量 的 情况 下 ,可 以 尝试 改进 对 
高 阶 窍 的 估计 方法 , 即 :寻找 一 个 适当 宽 一 点 的 Nch 
bin RRB A Nch, 这 样 每 个 Nch bin 中 的 事件 数 
就 增加 了 ,每 次 计算 xa? 的 精准 度 就 可 以 得 到 提高 。 


的 问题 是 :目前 和 未 来 实验 的 统计 量 是 否 能 保证 高 
阶 矩 测量 的 精准 度 , 或 者 说 多 少 统计 量 才能 保证 其 
精准 度 ? 


当然 ,适当 加 宽 Nch bin 的 条 件 必 须 保 证 不 影响 其 代 
表 的 初始 体积 , 即 :必须 与 每 个 Nch 有 相同 的 分 辨 
率 。 为 了 寻找 适当 的 Nch bin 宽度 ,将 0%~5% 中心 


另外 ,相对 论 重 离子 碰撞 中 ,守恒 荷 高 阶 矩 测量 
方法 的 特殊 性 也 使 我 们 必须 关注 其 统计 量 的 基本 要 
求 。 在 高 阶 矩 测量 值 的 估计 中 , 其 对 应 碰撞 中 心 度 
参数 由 多 重 数 区 间 诀 定 , 为 了 消除 这 个 区 间 中 不 同 
E 数 所 代表 初始 体积 的 涨 落 , 对 给 定 中 心 度 的 多 
重 数 区 间 , 人们 首先 在 每 个 多 重 数 Nch 计算 其 高 阶 
抢 ,然后 把 每 个 Nch 的 事件 数 做 为 权重 因子 ,对 应 中 
心 度 的 高 阶 矩 是 Nen X M) Cbin) P EA Neh ie BT FE 
的 加 权 平 均 ,这 方法 就 是 所 谓 的 中 心 度 bin 宽度 方 
YF, Bl: CBWC AVE. 

显然 ,这 方法 会 导致 一 个 问题 , 那 就 是 虽然 实验 
的 总 统计 量 都 是 百 万 量 级 的 事件 数 , 而 分 配 到 每 个 
Nch bin( 特 别 是 0%~5% 的 Nch bin) ,甚至 每 个 Nch 
上 的 统计 量 却 非 常 有 限 。 以 AutAu fb i Vsw= 
7.7 GeV 为 例 , 目 前 用 来 分 析 净 质子 xo? 的 总 事件 数 
(0%~80% 中心 度 ) 有 3 M。 在 0%~5% 中 心 度 , 确 定 
中 心 度 的 Neh bin 至少 有 100 个 Nch, 粗 略 平均 ,每 个 
Nch 上 只 有 1 900 个 事件 ,所 以 每 个 Ncn 上 计算 得 到 
的 高 阶 和 矩 的 统计 误差 非常 大 。 
因此 ,由 于 守恒 荷 高 阶 矩 本 身 的 特点 ,以 及 其 测 
量 方法 的 特殊 性 ,必须 估计 其 精确 测量 所 需要 的 基 
本 统计 量 。 下 面 先 以 四 阶 矩 为 例 进行 讨论 2。 我 
们 知道 ,只 有 在 统计 量 足 够 多 ,实验 随机 观察 量 的 分 
布 才 是 正 态 分 布 , 即 :满足 中 心 极 限定 理 CCentral 
Limit Theorem, CLT)。 只 有 在 满足 CLT 的 条 件 下 ， 
观测 量 的 平均 值 才 等 于 其 真实 值 , 并 且 不 依赖 于 统 
ihe, BCX), = CX), = A’ R P nnn NA 
同 的 统计 量 。 

我 们 用 UrQMD 模型 模拟 金 - 金 碰撞 ,在 Vsw= 
11.5 GeV 能 量 下 ,总 共产 生 了 250 M 个 事件 ,这 个 寻 
件数 的 量 级 与 目前 RHIC 相对 论 重 离子 碰撞 实验 的 
统计 量 相 当 。 然 后 ,将 总 样本 分 别 随 机 分 成 250、50 
及 10 个 子 样本 ,每 个 子 样本 对 应 的 统计 量 分 别 为 
1 M、5 M 及 25M。 在 每 个 子 样本 中 ,用 CBWC 方 法 


lirli 
Th 


O 
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度 分 成 不 同 的 bn 宽度 ,如 25 个 bin(60.2% 中 心 度 )、 
5 个 bin(61.0% Fly FE) 2 4 bin (62.5% 中心 度 )、 每 
个 Nch 以 及 直接 用 0%~5% 中心 度 。 同 时 ,为 了 讨论 
在 不 同 Nch bin 宽度 情况 下 所 得 结果 对 统计 量 的 依 
赖 性 ,我 们 把 模型 产生 的 数据 随机 分 成 250、125、 
50、25、10 以 及 5 个 子 样本 ,每 个 子 样本 对 应 的 统计 
量 分 别 为 :1 M、2M、5M、10 M、25M 以 及 50 M。 
按照 CBWC 方 法 ,首先 计算 每 个 子 样本 对 应 的 ka?， 
然后 平均 得 到 《xa?》, 对 于 净 质 子 和 净 电 蓓 的 结果 分 
别 如 图 1(a) 和 (Cb) 所 示 。 
由 图 1 可 见 , 当 总 统计 量 大 于 1M 后 ,用 51.0%、 
62.5%、65.0% 得 到 的 《xo?) 都 与 统计 量 无 关 , 但 是 ， 
对 于 每 个 Nch, 只 有 当 总 统计 量 大 于 25 M 以 后 ,所 
得 结果 才 与 统计 量 无 关 , 且 与 用 51.0% 的 结果 一 致 ， 
表明 此 时 51.0% 可 以 与 每 个 Nch 一 样 准 确 地 代表 初 
始 体 积 。 同 时 ,我 们 看 到 62.5% 及 55.0% 得 到 的 
《ko 明显 要 高 于 每 个 Nch( 或 51.0%) 计 算 的 结 
说 明 这 样 选择 的 Nch bin 宽度 太 大 了 ,不 能 与 每 个 
Neh 一 样 代 表 初 始 体 积 , 因 此 ,只 有 51.0% 可 以 代替 
每 个 Nch。 
可 见 , 对 于 四 阶 矩 ,基本 的 统计 量 要 求 是 :如 果 
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CBWC Method at 
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1 不 同 的 CBWC bin 宽度 情况 下 得 到 的 净 质 子 与 静电 荷 

的 对 xa* 统 计量 的 依赖 性 

Fig.1 Statistics dependence of xo’ of net-proton (a) and net- 

charge (b) multiplicity distributions calculated by the CBWC 
method at various centrality bin width 
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确定 QCD 相 边 界 的 若干 问题 


按 每 个 Nch 计算 ,总 统计 量 必须 大 于 25 M; 如果 按 
61.0% 计算 ,总统 计量 到 1 M 量 级 就 可 以 了 ,并且 ,该 
结论 其 他 碰撞 中 心 度 情况 下 也 成 立 。 因 此 ,在 总 样 
本 统计 量 低 于 25 M 的 情况 下 ,由 每 个 Nch 计算 得 到 
的 四 阶 矩 很 可 能 不 可 靠 , 利 用 51.0% 中 心 度 作 为 
Neh bin 宽度 计算 的 结果 会 更 可 靠 。 

对 于 净 质 子 的 六 阶 和 矩 *” 忠 ,基本 统计 量 的 要 求 
比 四 阶 矩 更 高 。 我 们 研究 了 在 0.4 GeV:c! 
0.8 GeV -c !' 以 及 |y|<0.5 截断 条 件 下 , 当 总 事件 统计 
量 从 1 M 增 加 到 300 M 时 ,在 0%~30% 中 心 度 ,每 个 
Neh 下 计算 得 到 的 (C6/C,) 一 直 随 着 统计 量 的 增加 
而 增加 ,说 明 在 这 个 相 空间 截断 下 ,即使 统计 量 达 到 
300 M 量 级 ,测量 得 到 的 CyC;, 都 不 可 靠 。 
1.2 ”统计 涨 落 的 影响 

我 们 知道 , 泊 松 分 布 可 以 描述 相互 独立 \ 彼 此 没 
有 关联 粒子 的 随机 涨 落 , 即 统计 涨 落 。 如 果 一 个 随 
机 变量 X, 满 足 泊 松 分 布 , 且 其 平均 值 为 1, 那么 ,对 
于 随机 变量 为 三 0,1,2,... 非 负 整 数 , 概 率 分 布 为 
f(k; A)=Pr(X=k) =e YE! ,RHIC 实验 中 产 
生 的 粒子 ,分 别 有 带 重子 数 为 1 的 强 子 ,带电 荷 数 为 
1 或 者 2 的 强 子 ,也 有 带 奇 异 数 为 1.2 或 3 的 强 子 ,以 
及 它们 相应 的 反 粒 子 。 从 最 简单 的 情形 出 发 ,假设 
重子 数 CX2 和 反 重 子 的 数 CV2) 都 服从 泊 松 分 布 。 
那么 PETAN, = NP - V2) 的 概率 分 布 是 两 个 
泊 松 分 布 的 互 关 联 (cross-correlation , E} : 


J(wi(NP)(N))- Df (Nat 
ETE 
(NENE S SNAN) OO 

Hh (NP) ACN?) OF Sl E NP R NA, 的 平均 值 ; 

Jy(z) 是 修正 后 的 第 一 类 贝 塞 尔 函 数 。 式 (2) 是 一 

个 标准 的 Skellam 分 布 汪 ,这 样 给 出 的 净重 子 数 概率 

分 布 与 HRG 模型 得 到 的 结果 相同 入 。 净 重子 的 高 

It KE Cac’) HY FH FE EY) AE A PRI 28 (Cumulant-generating 

Function, CGF ) 得 到 : 

Ki T (3) 

其 中 : Ks(t; (NEY, (NS))= InG(e'; (NP), 

(N2))sG(ts (NP), (N2,)) Æ Skellam 分布 的 概率 

PGF, 即 : 


<p < 


G(t; (N?),(N2,)) = Di lms (NEY, (NE Y) 
= G(t; (NEY )G(1/t; (N2,)) 


= eT (ONE) + (NE) + (nt et CE it (4) 

因此 ET E PAB PIN : 

Ke = (NI) +(N2,), ee = (NP) - (NA) (5) 

它们 由 重子 和 反 重 子 数 的 平均 值 唯一 确定 。 对 
于 带电 粒子 、 带 奇异 荷 的 粒子 ,如 果 每 种 粒子 的 多 重 
数 服从 泊 松 分 布 ,我 们 可 以 类 似 地 推导 ,它们 也 都 和 
玻 尔 效 曼 近似 下 的 强 子 共振 气体 模型 得 到 的 结 
二 和。 

根据 RHIC 三 个 碰撞 能 量 (Vsw=19.6 GeV. 
62.4 GeV 和 200 GeV) ,9 个 中 心 度 下 (能 量 为 
19.6 GeV 时 ,只 有 5 个 中 心 度 ) 的 质子 和 反 质子 数 的 
平均 值 ,我们 分 别 计算 了 净 质 子 数 统计 涨 落 的 标准 
差 o\ 偏 5S, 和 峰 度 x,, 结 果 发 现 :所 有 的 统计 涨 落 的 
和 抑 都 与 实验 测 到 的 矩 靠 得 非常 近 P ,它们 之 间 的 差 
直 比 原 值 小 了 一 个 数量 级 。 因 此 ,RHIC 上 产生 的 净 
质子 数 高 阶 矩 的 中 心 度 和 能 量 依赖 性 主要 是 统计 涨 
落 的 结果 ,观察 量 守恒 荷 高 阶 矩 中 的 统计 涨 落 不 可 
忽略 ! 
1.3 ”混合 事件 方法 


对 于 一 个 真实 的 样本 ,各 种 非 临 界 效应 相互 联 
系 ,很 难 单独 地 将 某 一 种 效应 扣除 ,更 好 的 扣除 思路 
是 用 包含 相同 背景 的 本 底 样 本 , 即 :混合 事件 样本 。 
通常 ,一 个 混合 事件 中 的 粒子 是 ae 
始 事 件 中 挑选 ,同时 保持 混合 事件 样本 的 总 粒子 数 
与 原始 样本 一 致 。 因 此 , 混合 事件 样本 具有 与 原始 
样本 相同 的 统计 涨 落 、 初 态 系统 尺度 涨 落 .中心 度 分 
饼 宽 度 、 探 测 效率 以 及 实验 截断 。 
相对 重 粒子 对 撞 中 ,混合 事件 方法 已 经 被 应 用 
于 多 种 观测 量 , 例 如 两 粒子 关联 、 粒 子 产 额 \ 横 动量 
谱 、 椭 圆 流 等 。 不 同 的 观测 量 构 造 混合 事件 的 方法 
并 不 完全 相同 。 例 如 ,两 粒子 快 度 关联 ,关注 的 动力 
学 特征 是 在 给 定 的 事件 中 粒子 的 快 度 方位 ,所 以 其 
混合 事件 是 将 一 个 给 定 事件 中 所 有 粒子 的 真实 快 度 
方位 进行 随机 替换 。 对 于 守恒 荷 的 高 阶 矩 ,关注 的 
动力 学 特征 是 守恒 荷 的 粒子 数 。 所 以 混合 事件 的 目 
的 是 消除 守恒 荷 粒子 之 间 的 关联 ,以 及 给 定 守恒 荷 
粒子 数 与 事件 之 间 的 关联 。 

我 们 发 现 最 好 的 方法 是 最 随机 性 的 池子 方法 ， 
如 果 原 始 样 本 的 统计 量 足 够 大 , 即 大 于 1 000 个 事件 
时 ,粒子 总 数 可 以 近似 为 无 限 多 ,可 以 将 全 部 粒子 放 
入 一 个 池子 里 面 ,然后 ,从 池子 里 随机 挑选 入 .个 带 荷 
粒子 ,这 N. 个 粒子 可 以 近似 看 作 是 从 不 同 的 原始 事 
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件 中 选 出 ,它们 之 间 的 关联 近似 为 零 。 从 池子 中 随 
机 挑选 一 个 带电 粒子 的 概率 是 一 个 常数 ,混合 样本 
中 的 守恒 荷 粒 子 数 的 平均 值 与 原始 样本 保持 一 致 。 
这 样 得 到 的 混合 样本 不 再 具有 与 真实 事件 相关 的 临 
界 涨 落 , 但 是 ,保留 了 事件 的 整体 与 系统 特征 。 如 果 
将 动力 学 累积 矩 定义 为 原始 样本 的 累积 矩 减 去 混合 
样本 的 累积 矩 , 那 么 ,动力 学 累积 矩 就 代表 了 临界 相 
关 的 涨 落 。 

利用 多 相 输 运 模 型 CA Multiphase Transport 
Model,AMPT) 的 default 模型 模拟 产生 了 19.6 GeV 
金 - 金 碰撞 样本 ,针对 统计 涨 落 、 中 心 度 分 饼 宽度 ,以 
及 探测 器 效率 的 扣除 公式 ,逐一 比较 动力 学 累积 
与 使 用 扣除 公式 后 的 结果 。 其 中 , 横 动量 截断 和 快 
度 截 断 分 别 为 :0.4GeV < p, < 0.8 GeV 和 -0.5< 
y < 0.5, 与 RHIC/STAR 实验 中 所 采用 的 截断 一 致 。 
质子 数 近似 为 总 净重 子 数 , 即 :N=N,。 

图 2 展示 了 xa** 的 中 心 度 依赖 性 ,为 了 观察 探 
测 器 效率 的 影响 ,分 别 以 80% 和 60% 的 概率 来 随机 
选择 粒子 ,然后 分 别 计算 对 应 的 累积 矩 。 其 中 圆圈 
标示 100% , 方 框 标示 80%, 三 角形 标示 60%。 图 2 
Ca) 是 在 三 种 不 同 的 探测 器 效率 影响 下 ,原始 样本 中 
xo” 的 中 心 度 依赖 性 ,可 见 , 三 组 点 之 间 相 互 分 离 ， 
累积 矩 随 探 测 器 的 效率 而 变化 , 且 探 测 器 效率 的 影 
响 依 赖 于 中 心 度 。 图 2(b) 是 在 三 种 不 同 探测 器 效 
率 影 响 下 , 混合 样本 的 xcz2" 的 中 心 度 依赖 性 ,探测 
器 效率 同样 影响 了 混合 样本 的 抢 ,而 且 与 图 2(Ca) 中 
对 于 原始 样本 的 影响 类 似 。 

利用 文献 [45-46j] 中 的 公式 ,图 2(d) 给 出 了 公 
式 修正 后 不 依赖 于 探测 器 效率 的 xo?2。 由 图 2(d) 可 
JL , 方 框 和 三 角 都 与 对 应 的 圆圈 在 误差 棒 的 范围 内 
重合 ,这 说 明 公式 修正 后 的 累积 矩 很 好 地 扣除 了 探 
测 器 效率 的 影响 。 图 2(c) 是 在 三 种 不 同 探测 器 效率 
影响 下 ,动力 学 ko 如 随 中 心 度 的 变化 ,可 见 三 组 不 同 
的 点 在 误差 范围 内 都 基本 重合 。 这 表明 动力 学 kai 
与 公式 修正 后 的 结果 一 样 ,都 不 依赖 于 探测 器 效率 ， 
都 扣除 了 探测 器 效率 的 影响 。 而 且 ,动力 学 累积 
不 仅 减 去 了 统计 涨 落 ,同时 还 减 去 了 中 心 度 和 探测 
器 效率 的 影响 。 


2 ”有 限 系 统 尺 度 的 影响 


2.1 ”有 限 系统 尺度 对 相 变 涨 落 的 影响 


在 大 重子 化 学 势 以 及 手 征 极 限 附近 ,由 于 格 点 
QCD 的 计算 还 存在 重重 困难 ,为 了 计算 有 限 尺 度 系 
统 的 热力 学 量 在 相 边 界 附近 的 涨 落 , 人 们 常常 借助 
一 些 唯 象 模型 和 有 效 理论 ,与 QCD 相 变 同属 一 个 普 
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2 探测 效率 分 别 为 100%、80% 和 60% 时 ,由 原始 样本 
(a)`、 混 合 样 本 (b)、 动 力学 (c) 和 公式 修正 后 (d) 所 得 到 的 xa? 随 
中 心 度 的 变化 

Fig.2 The centrality dependency of xo” of original (a) and 
mixed (b) samples, dynamical (c) and formula corrected (d) for 

three values of efficiency 100%, 80% and 60% 


TR AY BEL = ERA OAM OO) H JERE , DI 
及 三 维 三 态 Potts 模 型 。 三 维 Ising 模 型 是 由 N 个 排 
列 在 LxLxXL 个 晶 格 点 阵 上 的 粒子 组 成 ,粒子 的 自 旋 
方向 可 以 向 上 ,也 可 以 向 下 ,用 +1 表 示 向 上 ,-1 表示 
向 下 ,系统 的 状态 可 以 由 一 系列 自 旋 表示 , 即 : 
{51,52， Sps; = E1, i= 1,2,7, N, A FIE f 
符 {s} 表 示 。 如 果 仅 考虑 最 近邻 粒子 的 相互 作用 ， 
以 及 粒子 与 外 部 磁场 五 的 相互 作用 ,总 能 量 可 以 表 
AA: 


N 
E) =J Xss; = H Y's: (6) 


其 中 :符号 < ij > 表示 只 对 所 有 最 近邻 粒子 的 自 旋 
求 和 。 系 统 在 某 个 微观 状态 { s,} 的 平均 总 磁化 强度 
和 平均 磁化 强度 分 别 为 : 


当 温 度 太 7<7, 自 旋 - 自 旋 相 互 作用 相对 较 强 ,所 
有 的 粒子 倾向 于 呈 同 一 指向 排列 ,此 时 的 | 四 | 接近 于 
1, RATATA; SBE k TJ, 自 旋 - 自 旋 相 互 作 
用 相对 较 弱 ,每 个 格 点 的 自 旋 状 态 随机 地 向 上 ,或 者 
向 下 ,|m| 接 近 零 ,系统 处 于 完全 无 序 的 状态 。 在 没 
有 外 场 的 情况 下 ,系统 从 高 温 无 序 相 到 低温 有 序 相 
的 转变 是 连续 相 变 , 相 变 点 称 为 居 里 点 , 居 里 点 是 温 
度 T 和 外 磁场 且 平 面 上 的 临界 点 ,临界 点 的 温度 用 
表示 。 人 们 通过 有 限 尺度 标 度 分 析 方 法 ,计算 出 
了 三 维 伊 辛 模型 的 临界 温度 7=4.51” ,与 重 整 化 群 
理论 的 结果 一 致 ”。 
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三 维 三 态 Potts 模 型 是 有 限 温 度 纯 规范 QCD 理 
论 的 一 个 简单 有 效 模型 , 它 的 临界 点 、 相 边界 可 以 通 
过 三 维 Ising 模 型 直接 映射 到 QCD 相 边 界 。 因 此 ， 
我 们 利用 三 维 三 态 Potts 模 型 ,以 及 其 磁化 强度 高 阶 
和 矩 行为 的 普 适 性 ,展示 其 磁化 强度 的 高 阶 感应 率 在 
一 级 相 变 、 临 界 点 、 平 滑 过 渡 区 域 的 行为 ,以 及 这 些 
高 阶 感应 率 的 有 限 尺度 特征 。 
Potts 模型 每 个 格 点 上 的 自 旋 可 以 取 多 个 离散 
值 ,通常 这 些 值 用 从 1 开始 的 整数 表示 ,一 个 g 态 的 
Potts 模 型 ,其 自 旋 可 以 取 整 数 s=1,2,… ,gq。 任 意 相 
邻 的 两 个 自 旋 ,如 果 有 相同 的 值 , 则 贡献 一 个 7 给 哈 
密 顿 量 ;如 果 不 相同 , 则 贡献 0。 如 果 自 旋 刚 好 和 外 
场 的 方向 相同 , 则 贡献 及 给 哈密 顿 量 , 如果 不 相同 ， 
则 贡献 0。 因 此 哈密 顿 量 可 以 表达 如 下 : 
H’ =- JE ðh s» 8) - HEd(s,, S) (8) 
Hp JERI H DE (ay) ZT A TAE FA, ER 
们 让 其 值 为 1。 符 号 5(s,,s,) 在 i=j 时 取 1,iz#j 时 取 
0。s, 是 所 谓 的 ghost spin, 它 的 方向 就 是 外 场 的 方 
向 Potts 模型 的 配 分 函数 可 定义 为 Z(p,h)= 
Xeno, ERE = 3.6(s,,s,) 是 相 邻 自 旋 之 间 
{s,} (i) 


的 相互 作用 能 


;M = >5(s,s。) 是 磁化 强度 ;有 = BH 


是 约 化 外 场 ;8 = 1/T 是 温度 的 倒数 。 
三 维 三 态 Potts 模 型 的 自 旋 s, 可 以 取 1、2、3 三 种 
状态 ,并 且 自 旋 位 于 一 个 立方 唱 格 上 。 满 足 Z(3) 对 


称 性 ,Z(3) 对 称 性 的 自发 破 缺 被 认为 是 导致 济 火 
QCD 的 退 禁 闭 相 变 的 原因 ,而 加 入 外 场 导致 2(3) 对 
称 性 破 缺 就 像 QCD 里 面 费 米子 多 项 式 破坏 了 Z(3) 


中 确定 QCD 相 边 界 的 若干 问题 
型 的 序 参量 可 用 磁化 强度 的 平均 值 = 3p f 
来 定义 51。 


在 QCD 纯 规 范 理论 里 , 序 参量 是 Polyakov 
loop , 它 起 到 了 磁化 强度 的 作用 ,表征 退 禁 闭 相 的 中 
心 对 称 性 自发 破 缺 。 三 维 三 态 Potts 模型 在 零 外 场 
的 情况 下 是 一 个 温度 驱使 的 一 级 相 变 ， oo 
渐 增 大 时 ,一 级 相 变 将 变 弱 并 终止 于 一 个 二 级 相 变 
点 , 即 临 界 点 ,越过 临界 站 
这 个 与 QCD 相 变 类 似 ,QCD 中 重子 化 学 势 的 作用 
相当 于 Potts 模型 中 外 场 的 作用 ,它们 的 变化 导致 温 
度 驱 使 的 相 变 从 一 级 变 到 二 级 ,再 到 平滑 过 渡 。 

对 于 一 级 相 变 、 二 级 相 变 和 平滑 过 渡 , 分 别 选 取 
了 外 场 h=0.000 5、0.000 775、0.000 2 三 个 值 。 令 对 = 
《M), 在 每 个 外 场 下 ,模拟 了 尺度 LL 在 40~70 之 间 不 
同 大 小 的 系统 ,分 别 计算 了 二 阶 感应 率 (x, = C/V), 
给 出 了 它们 在 相 变 温度 附近 的 行为 。 三 个 外 场 下 的 
二 阶 感应 率 分 别 在 图 3 的 三 个 子 图 中 显示 出 来 。 每 
个 子 图 中 不 同 颜色 的 线 代 表 不 同 的 系统 尺度 。 从 图 
中 可 以 看 出 ,在 一 级 相 变 (图 3(a))、 二 级 相 变 (图 3 
(b)) 以 及 平滑 过 渡 ( 图 3(c)) 的 时 候 , 二 阶 感应 率 在 
相 变 温度 (7,) 附 近 都 呈现 出 一 个 峰 。 我 们 知道 ,在 
热力 学 极限 下 ,对 于 一 级 相 变 ,二 阶 感应 率 应 该 是 一 
个 5 函数 ,对 于 二 级 相 变 ,二 阶 感应 率 会 发 散 。 有 限 
的 系统 使 它 在 两 种 相 变 中 都 变 成 一 个 有 限 的 
峰 c25。 随 着 系统 尺度 的 增加 , 峰 的 位 置 向 高 温 方 
向 移动 。 比 较 图 3 中 的 三 个 子 图 ,可 以 看 出 ,在 相 变 
温度 附近 ,不 同 的 外 场 下 ,二 阶 感应 率 的 基本 结构 是 
相似 ,都 有 一 个 峰 。 表 明 在 相 变 温度 附近 ,有 限 系 统 
感应 率 会 出 现 非 单调 行为 ,这 些 非 单调 行为 不 仅 发 


对 称 性 一 样 。 文 献 [49-50j 指 出 ,在 高 温 重 夸克 质量 。” 生 在 二 级 相 变 的 临界 温度 附近 ,也 会 发 生 在 一 级 相 
的 情况 下 ,QCD 可 简化 为 三 维 三 态 Potts 模 型 。 该 模 。 ” 变 和 平滑 过 渡 区 域 。 
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图 3 在 4 个 不 同 尺度 L=40, 50, 60, 70 的 系统 中 ,不 同 的 外 场 h=0.000 5 (a), 0.000 775 (b), 0.002 (c) F , 


二 阶 感应 率 x, 随 温度 7/7, 变 化 的 行为 
Fig.3 x, vs. T/T „at h= 0.000 5 (a), 0.000 775 (b), 0.002 (c) with the system sizes of L = 40, 50, 60, 70 


2.2 ”有 限 系 统 临 界 参 数 的 确定 


相 图 给 出 了 发 生 相 变 边 界 的 位 置 , 为 了 精确 确 
定 相 边界 的 位 置 ,必须 先 给 出 一 些 关 键 点 的 位 置 ， 


如 :二 级 相 变 的 临界 点 就 是 一 个 这 样 的 关键 点 ,临界 
点 的 温度 到, 即 临界 温度 就 是 人 们 最 感 兴趣 的 参数 。 
我 们 知道 ,如 果 系 统 的 体积 足够 大 ,可 由 磁化 率 分 
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布 的 峰值 位 置 近似 确定 ;但 是 ,从 上 面 的 讨论 我 们 知 
道 ,对 于 有 限 系 尺度 系统 ,其 峰值 的 位 置 随 系 统 尺度 
变化 。 如 果 通 过 有 限 系 统 观 测量 的 涨 落 来 确定 相 变 
参数 ,那么 ,必须 利用 有 限 尺度 标 度 律 。 

对 于 有 限 尺度 系统 ,如 果 发 生 二 级 或 一 级 相 变 ， 
ne ee eh y 


= 


结果 可 以 标 度 到 一 个 普 适 的 等 式 , 这 就 是 所 谓 的 有 
限 尺 度 标 度 律 , 即 : 
Q(T, L) =L” foi") (9) 


SUH: OCT, L) Fe WM E , — A A i EA RARE 
PRA; = CT-T,./T, AZ Hoda RE , T, HAL BT AT BE 
到 的 相 变 温度 :为 标 度 函 数 ; ;4 和 v 分 别 是 观察 量 O 
和 关联 长 < 对 应 的 标 度 指数 。 其 中 
XWv, 对 于 二 级 相 变 临界 点 ,a 是 介 于 4d~0 之 间 的 分 数 ; 
对 于 一 级 相 变 ,a 只 依赖 于 系统 维度 d, 是 维度 的 整 
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数 倍 (n-1)q, 其 中 代表 自由 能 对 温度 和 外 场 偏 导 
的 阶 数 。 

图 4(a) 给 出 了 在 无 外 场 的 情况 下 ,二 维 Ising 模 
型 二 级 相 变 的 有 限 尺度 标 度 行为 ,图 4(c) 是 在 外 场 
h=0.000 5 的 情况 下 ,三 维 三 态 Potts 模 型 一 级 相 变 的 
有 限 尺度 标 度 行为 。 由 图 4 可 见 , 不 同 尺度 标 度 的 
磁化 强度 曲线 在 临界 点 附近 和 一 级 相 变 区 域 都 很 好 
地 重合 在 一 起 。 通 常 ,人 们 利用 这 个 标 度 规律 确定 
临界 温度 和 临界 指数 时 ,都 是 预先 给 定 这 两 个 参数 ， 


然后 看 标 度 曲线 是 否 能 重 琶 在 一 起 , 即 :通过 观察 不 
同 尺度 标 度 的 观测 量 曲 线 在 约 化 温度 处 的 重合 情 


况 ,评估 事先 假定 的 相 变 参数 是 否 正确 ”*。 显 
然 , 这 种 方法 存在 很 大 的 不 确定 性 ， 标 度 曲 线 
重合 情况 大 多 由 肉眼 判断 ,无 法 定量 描述 ;另外 ,在 
相 变 约 化 温度 附 近 , 相 应 的 重合 区 域 比较 天 ,如 图 4 
Ca) 和 (c) 所 示 , 加 之 标 度 尺 度 因 子 的 加 入 ,一般 无 法 
精确 确定 对 应 的 相 变 参数 。 
如 果 将 标 度 曲线 的 自 变量 选 为 温度 ,那么 , 标 度 
曲线 在 偏离 相 变 温度 时 就 会 彼此 分 开 , 如 图 4(b) 和 
Cd) 所 示 , 仅 在 相 变 温度 处 , 标 度 曲线 交 于 一 点 , 即 固 
定点 。 固 定点 对 应 的 温度 和 标 度 指数 率 即 为 相 变 温 
度 和 相 变 标 度 指数 率 。 显 然 , 相 对 于 双 标 度 曲 线 ( 图 
4(a) 和 (c)), 单 标 度 曲线 的 固定 点 (图 4(b) 和 (d)) 能 
够 完全 精确 地 给 出 对 应 的 相 变 参数 , 而且, 定量 判 定 
国定 点 比 曲 线 重合 要 简单 、 容 易 得 多 。 
根据 重 整 化 群 理论 ,临界 点 为 不 稳定 固定 
点 555, 一 级 相 变 线 上 是 不 连续 的 固定 点 震中 。 在 平 
滑 过 渡 相 转变 a 于 系统 尺度 ,如 果 


将 有 限 尺 度 标 度 性 推广 到 平滑 过 渡 区 域 ,这 时 标 度 
"o aA 
图 4(b) 和 4(d) 可 以 看 到 , 当 温 度 越 来 越 远 离 


图 4 ” 标 度 平均 磁化 强度 mL” 随 自 变 量 标 度 约 化 温度 碟 * 
(a, c) 和 温度 7(b, d) 的 变化 关系 

(a, b) 无 外 场 h 且 0 下 的 二 维 Ising 模 型 ,(c, d) 外 场 h=0.000 5 下 

的 三 维 三 态 potts 模 型 

Fig.4 The scaled order parameter mL” vs. the scaled variable 

tL (a, c) and T (b, d) 
(a, b) 2D Ising model with external field h=0, (c, d) 3D three- 
state Potts model with h=0.000 5 

临界 温度 时 ,不 同 尺度 标 度 的 观测 量 也 越 来 越 分 散 ， 

反之 , 当 温度 越 来 越 靠 近 临 界 温 度 时 ,不同 尺度 标 度 

的 观测 量 逐 渐 聚 拢 ,在 临界 温度 处 汇聚 成 固定 点 。 

因此 ,可 以 根据 不 同 尺 度 观 测量 的 聚集 度 ( 分 布 的 宽 

度 ) 来 定义 固定 点 。 在 固定 温度 和 标 度 指数 率 时 , 定 

义 分 布 宽 度 为 不 同 尺 度 标 度 观测 点 OQ (T, L) LEAS 

的 标准 方差 的 平方 根 “, 即 : 
Asouor 

D(T,a)= Wee (10) 

式 中 :D(7, a)ren TREE BCE a HRN, 7E 

不 同系 统 尺度 的 个 数 ; AS ,是 所 有 标 度 观测 量 
相对 于 它们 期 望 值 的 偏差 加 权 方 差 值 , 即 : 

2 
N, T,L,) L} = T,L)L* 
2 区 Ae =O LII] 


O(TL)L 


(11) 


权重 因子 w,= 5[2(7 L) L] ÆR REPR E KAN 
Re Q(T, L) Li WR. (OCT, L) L) J D H 
望 值 。 

当 了 和 a 越 来 越 远离 相 变 值 时 ,不 同 尺度 的 标 度 
曲线 越 来 越 分 散 ,相应 D(7, a) 增 加 ; 当 7T 和 a 逐渐 
接近 相 变 值 时 ,不同 尺度 的 标 度 曲线 逐渐 聚拢 ,相应 
D(T, a) 减 小 。 当 7 和 a 都 达到 相 变 值 时 ,理论 上 不 
同 尺度 的 标 度 点 0O(7,L)L 取 相同 值 ,并 相交 于 固 
定点 ,这 时 D(7T, a) 达 到 极 小 值 D(7, a)， 约 为 1。 
随 着 温度 和 标 度 指 数 率 的 改变 , 标 度 观 测 点 的 分 布 
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宽度 呈现 分 散 和 聚拢 的 趋势 是 形成 固定 点 的 必要 条 ” 上 , 即 :T 和 a Rg de cli 
件 。 由 于 观测 量 的 测量 误差 ,以 及 其 他 实际 测量 过 ” ”的 分 布 宽度 D(7, a), 如 图 5 所 示 的 等 高 线 。 图 5 每 
程 中 的 不 确定 因素 ,固定 点 并 不 是 一 个 理想 的 点 , 实 。 个 子 图 右 侧 的 颜色 线条 代表 宽度 D(7T, a) 的 大 小 ， 
际 上 相应 D(7T, a), 会 比 1 大 。 其 中 :红色 和 蓝 色 分 别 代 表 最 小 值 和 最 大 值 。 在 了 
如 果 改 变 T 和 a,D(7, a) 存 在 极 小 值 ,并 且 约 为 ” 和 a 平面 内 ,各 类 颜色 分 布 代表 标 度 曲 线 分 布 宽度 
1, 就 意味 着 标 度 观察 量 O(7,L)L 曲线 相交 于 固定 ”也 在 不 同位 置 的 大 小 。a 的 取 值 范围 为 -1.2~1.15 ,其 
点 ,固定 点 所 对 应 的 T 和 a 就 是 相 边界 的 特征 参数 。 ”中 间隔 为 0.05。 
不 同 相 变 类 型 对 应 的 特征 参数 不 同 ,临界 点 处 标 度 当 外 场 取 /=0.000775 和 0.0005 时 ,对 应 图 5 
指数 率 为 分 数 ,一 级 相 变 线 上 标 度 指 数 率 是 与 系统 (a) 中 的 三 级 相 变 点 和 图 5(b) 一 级 相 变 线 附近 。 通 
维度 有 关 的 整数 ,crossover 区域 由 于 观察 量 不 依赖 过 改变 T 和 a, 图 5(a) 和 (b) 中 的 力 值 会 出 现 逐 渐 减 
于 系统 尺度 ,在 整个 相 转 变温 度 区 域 D(7T, a) 都 近 ”小 的 一 块 红色 区 域 ,其 中 :颜色 最 红 的 位 置 束 是 宽度 
似 为 常数 1, 对 应 的 标 度 指数 率 a 为 0。 取 极 小 值 的 位 置 ,不 同 尺 度 标 度 的 序 参量 曲线 在 此 
现在 ,以 三 维 三 态 Potts 模 型 在 三 个 相 变 区 域 产 。” 交 于 一 点 , 即 : 固 定点 。 图 5(a) (或 (b)) 中 两 条 相互 
生 的 数据 为 例 ,说 明 如 何 应 用 上 面 对 固 定点 定量 描 ”垂直 点 线 的 交点 标 出 了 相 变 参数 平面 内 宽度 取 极 小 
述 的 方法 来 确定 相 边界 参数 。 三 个 相 变 区 域 分 别 是 ” 值 的 位 置 ,该 位 置 对 应 的 相 变 参数 束 是 临界 温度 和 
二 级 相 变 的 临界 点 .一 级 相 变 线 和 平滑 过 渡 区 ,对 应 标 度 指数 率 ( 或 一 级 相 变 的 温度 和 标 度 指 数 律 )。 因 
的 外 场 为 :=0.000 775 0.000 5 和 0.002 ,观测 量 为 平 。 此 , 相 变 参数 平面 内 的 宽度 等 高 图 精确 给 出 了 临界 
均 磁 化 强度 , 即 : 序 参量 。 首 先 ,在 相 变 参数 的 平面 。 区 域 和 一 级 相 变 固定 点 的 位 置 。 
(a) h= 0.000 775 Den Gas (b) h= 0.000 5 , r (c) h = 0.002 
70 0.50 


oo 


+} lilo , 0 06 
1.8195 1.8200 1.8205 1.8210 18180 18185 1.8190 1.8195 1.825 4 1.825 8 1.8262 
T T T 


图 5 三 个 外 场 处 ,在 温度 7 和 标 度 指数 率 a 的 平面 上 分 布 宽度 D(7,a) 的 等 高 图 
(a) h = 0.000 775 ,(b) A = 0.000 5, (c) A = 0.002 
Fig.5 Contour plot of D (T, a) for three different external fields (a) A = 0.000 775, (b) A = 0.000 5, (c) A = 0.002 

当 外 场 取 本 0.002 时 ,对 应 图 5(c) 中 的 平滑 过 导 的 能 量 样 的 观测 量 , 如 内 能 、 比 热 ,在 临界 ( 相 变 ) 
渡 区 域 ,这 里 力 值 的 等 高 线 平行 于 7T 轴 , 呈 带 状 分 ”区 域 受 非 奇 异 部 分 的 影响 ,其 有 限 尺度 标 度 行为 会 
布 ,带子 的 颜色 逐渐 变化 ,在 黑色 的 点 线 a=0 附 近 最 ”破坏 ,没有 明显 的 固定 点 行为 。 而 且 , 对 于 像 系 统 尺 
红 ,D(7 aq) 接近 常数 1, 这 些 特征 完全 不 同 于 图 5 ，” 度 、 温 度 样 的 物理 量 如 何 与 相对 论 重 离子 碰撞 中 的 
Ca) 和 (b) 中 的 分 布 。 表 明 在 平滑 过 渡 区 域 , 标 度 曲 “ 相关 物理 量 对 应 ,还 需要 认真 地 研究 决定 。 
线 不 依赖 于 系统 尺度 ,D(7, a) 不 依赖 相 过 渡 温 度 ， 


它 的 值 只 随 着 a 改变 , 相 过 渡 的 标 度 指数 率 为 0。 3 ”有 限 演化 时 间 的 影响 
由 此 可 见 , 分 布 宽度 值 D 的 大 小 对 相 边 界 上 的 ee 
固定 点 行为 十 分 敏感 ,可 以 很 好 地 描述 固定 点 行为 ， 通常 ,基于 格 点 QCD 和 各 种 QCD 模型 的 计算 


确定 相 边界 相 变 特征 参数 ,以 及 区 分 不 同 的 相 变 类 ”都 假设 系统 达到 了 热力 学 的 平衡 态 ””” 。 但 是 , 相 
型 。 当 然 , 具体 将 固定 点 方法 应 用 到 重 离子 碰撞 实 ”对 论 重 离子 碰撞 中 形成 的 火球 存在 的 时 间 非 常 有 
验 RHIC 能 量 扫描 计划 中 ,不 仅 需 对 诸如 观测 量 的 限 , 只 有 约 10 fm'c ,系统 很 可 能 没有 时 间 演 化 到 热 
选取 做 出 谨慎 的 讨论 ,因为 ,有 限 尺度 系统 的 标 度 规 pean ”。 一 般 可 能 的 状态 是 以 下 三 种 情 
律 和 固定 点 行为 只 对 某 些 像 自 由 能 对 外 场 的 偏 导 ， ” 况 :1) 如 果 碰 撞 能 量 足 够 高 ,系统 经 过 热 化 达到 解禁 
即 :磁化 强度 、 磁 化 率 , 会 表现 出 较 好 的 有 限 尺 度 标 。 闭 相 变 所 需 的 温度 , 系统 达到 热平衡 ,夸克 和 腰子 
度 性 和 明显 的 固定 点 行为 ;对 于 像 自由 能 对 温度 偏 等 离子 体 (Quark-Gluon Plasma,QGP) 形 成 ;2) 也 可 
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EE 系统 达到 局 域 热平衡 ,形成 QGP 的 droplet;3) 还 
[能 完全 没有 达到 热平衡 ,没有 经 历任 何 相 变 。 因 
L, 平 衡 \ 局 域 平衡 和 非 平衡 是 相对 论 重 离子 碰撞 后 
JT 能 形成 的 三 种 状态 。 此 外 ,在 临界 区 域 ,临界 慢 化 
也 会 使 系统 更 难于 达到 热力 学 平衡 态 。 

通常 人 们 利用 动力 学 演化 方程 (例如 Langevin 
B) JIEN FESO 流体 力学 方程 汪 % 和 各 种 的 弛 豫 模 
Foe 0) 来 估计 QCD 临界 区 的 非 平 衡 效应 。 但 
是 ,目前 动力 学 演化 方程 仅 在 临界 点 .平滑 过 渡 区 域 
能 精确 求解 ,在 一 级 相 边界 处 尚 无 法 精确 求解 。 我 
们 知道 ,Metropolis 算法 能 够 很 好 地 模拟 Ising 模型 
的 非 平 衡 演 化 ”*。 这 样 可 以 方便 地 得 到 整个 相 边 
界 上 非 平衡 弛 豫 过 程 的 数值 结果 。 然 后 ,给 出 Ising 
模型 的 非 平 衡 演 化 特征 ,并 引入 了 描述 非 平 衡 演化 
的 时 间 尺 度 。 最 后 ,讨论 非 平 衡 演化 对 观测 量 的 
影响 。 
3.1 ” 非 平衡 演化 的 时 间 描 述 和 特征 

从 一 种 初始 位 形 开始 ,Ising 系统 可 以 自发 演化 
到 平衡 态 ,达到 平衡 之 前 的 演化 称 为 弛 豫 过 程 或 非 


ZI 


= | 


z| 


平衡 演化 。 图 6(Ca~d) 给 出 了 从 随机 位 形 开始 ,系统 


. TIT,= 0.93 
(a) 
-0.2 


= 


尺度 L=60,4 个 温度 下 , 即 T/T. 43 Hl 0.93 0.99, 1.00 
和 1.03, 人 磁化 强度 的 绝对 值 |m| 随 时 间 的 演化 (使 用 m 
的 绝对 值 是 因为 在 h=0 时 磁化 强度 的 符号 随机 ,两 
条 曲线 是 从 样本 中 随机 选择 的 两 个 演化 过 程 ) ,每 个 
温度 进行 了 5 000 次 演化 。 前 两 个 温度 在 一 阶 相 变 
线 上 ,7/7T=1.00 是 临界 点 ,7/7T=1.03 是 平滑 过 渡 。 水 
平 轴 是 时 间 1, 即 ;sweep 数 (每 完成 一 次 对 所 有 格 点 
(N=LxLxL) 的 判别 , 称 为 一 次 扫描 (sweep) ,每 完成 
一 个 sweep ,系统 的 位 形 就 更 新 一 次 )。 两 条 曲线 是 
从 样本 中 随机 选择 的 两 个 演化 过 程 。 

在 图 6(a) 中 ,| 四 开始 逐渐 增加 ,然后 趋 于 稳定 
值 人 4/。 如 果 观 察 量 与 平衡 期 望 值 的 差异 大 于 平衡 时 
的 方 均 根 a = (x - 1)*)” ,系统 远 未 达到 平衡 。 从 
非 平衡 演化 到 平衡 的 时 间 称 为 弛 豫 时 间 。 在 模拟 
中 ,第 ;个 演化 过 程 的 弛 豫 时 间 世 是 在 四 | 值 进 入 区 
间 (4- ww+a) 后 的 步 数 。 网 6(a) 中 ,两 条 曲线 的 
稳定 值 相同 ,但 是 ,一 条 曲线 的 弛 了 豫 时 间 比 另 一 条 曲 
线 长 很 多 。 在 低温 时 ,不 同 曲线 的 弛 豫 时 间 差 异 显 
著 , 如 图 6(a) 和 (b); 在 高 温 时 ,差异 明显 减 小 ,如 图 
6(c) 和 (d) 所 示 。 
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6 温度 T/T 分 别 为 0.93 (a).0.99 (b).1.00 (c) FH 1.03 (d) 时 ,随机 初始 位 形 在 L=60 时 ,|m| 随 时 间 的 演化 过 程 
Fig.6 The evolution of |m| with time at L = 60 for random initial configuration at temperatures T/T, = 0.93 (a), 0.99 (b), 1.00 (c) and 
1.03 (d), respectively 
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式 中 :nn 是 演化 过 程 的 总 数 ;zi 是 第 i 个 演化 过 程 的 
弛 驳 时 间 。 显 然 ,在 zt, 时 ,并 不 是 所 有 系统 都 能 
到 平衡 态 , 其 中 部 分 系统 很 可 能 还 处 于 非 平 衡 状态 。 
通常 , 弛 豫 时 间 取 决 于 系统 演化 的 动力 学 机 制 、 
系统 尺度 ` 温度, 以 及 初始 位 形 等 。 图 7 给 出 了 平均 
att PR ANY E pE Ui RE. 系统 尺 度 和 初始 位 形 的 变化 。 在 
TBA EL ,平均 弛 锣 时 间 的 峰值 随 着 系统 尺度 而 增 大 ， 
如 图 7(a) 所 示 。 图 7(b) 表 明 ,无 论 选择 随机 初始 位 
形 还 是 极 化 初始 位 形 , 平 均 弛 驳 时 间 在 处 都 出 现 
峰值 ,这 是 由 于 临界 慢 化 所 导致 的 弛 豫 时 间 变 长 。 
理论 上 ,对 于 无 限 大 系统 ,在 临界 温度 附近 ， 
动力 学 的 弛 鸳 时 间 发 散 , 即 : 
roc oa |T-T,| (13) 


Ti (12) 


— 


式 中 :是 关联 长 度 的 临界 指数 ;z 是 动力 学 指数 , 它 
与 动力 学 的 普 适 类 相关 。 对 于 有 限 系统 尺度 , 弛 和 豫 
时 间 具 有 索 律 行为 , 即 : 
Tom cc 大 (14) 
Mp He, RATEI S PTE AP EY ts RANT Ta] Tog 
在 临界 温度 处 与 系统 尺度 的 双 对 数 图 ,其 线性 拟 合 
斜率 为 2.06+0.03 , 与 文献 [69-70] 一 致 。 因 此 ,平均 
弛 驳 时 间 云 也 随 着 系统 尺度 的 z 次 圭 发 散 , 与 ro 的 
行为 相同 。 这 表明 :模拟 系统 的 平均 弛 豫 时 间 与 动 
力学 弛 豫 时 间 一 样 , 也 能 够 很 好 地 量化 给 定 系 统 的 
时 间 演 化 特征 。 

对 于 7 和 ,平均 弛 豫 时 间 基 本 不 依赖 于 系统 尺 
度 和 初始 位 形 , 如 图 7(a) 和 (b) 所 示 , 无 论 系统 大 小 
和 初始 位 形 如 何 , 平 均 弛 和 鸳 时 间 都 很 短 ,并 趋 近 于 
零 。 这 是 由 于 接受 概率 4(u 一 v)( 式 (12)) 是 7 的 
增加 函数 ,温度 越 高 ,接受 概率 越 高 ,所 需 弛 驳 时 间 
越 短 ,系统 更 容易 从 一 种 位 形 演化 到 另 一 种 位 形 。 
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Random Config. 


~L=60 
+L=50 


L=60 


+ Random 
Polarized 


TIT, 


7 PIR E pE FZ AR A 
(a) 随机 初始 位 形 系统 尺度 L=50 和 60,(b) 固定 系统 尺度 L=60, 随 机 初始 位 形 和 极 化 初始 位 形 
Fig.7 Average relaxation time as a function of temperature 
(a) System sizes L = 50 and 60 starting from random initial configuration, (b) Random initial configurations and polarized initial 
configurations at a fixed system size L=60 


相反 ,对 于 T<7T ,平均 弛 耶 时 间 强 烈 依赖 于 系统 ”也 可 能 不 同 ,如 图 6(a) 中 的 曲线 所 示 。 所 以 ,如 果 测 
尺度 和 初始 位 形 。 首 先 , 图 7(a) 表 明 RARER ” 量 的 时 间 不 是 足够 长 ,那么 ,有 些 系 统 很 可 能 尚未 达 
大 ,平均 弛 殉 时 间 越 长 。 另 外 ,图 7(b) 还 表明 ,随机 ”到 平衡 ,观测 到 的 所 有 过 程 中 就 会 存在 一 定 比 例 的 
位 形 的 平均 弛 鸳 时 间 比 极 化 位 形 长 得 多 ,这 是 因为 非 平 衡 过 程 , 这 与 相对 论 重 离子 碰撞 的 情况 非常 
7< 工 的 平衡 态 接近 有 序 相 ,对 于 极 化 位 形 , 其 本 身 就 “相似 。 

是 有 序 相 ,因此 ,无 需 太 长 时 间 的 演化 就 可 以 达到 有 如 前 所 述 :相对 论 重 离子 碰撞 中 临界 涨 落 敏感 
序 相 的 平衡 态 ; 但 是 ,对 于 随机 位 形 ,其 本 身 接近 无 ”的 观测 量 是 序 参量 高 阶 矩 ,对 Ising 模型 来 说 , 令 式 
序 相 ,因此 ,需要 很 长 时 间 的 演化 才能 达到 有 序 相 的 。 CL AX = ||, 6X = 四 | - 《| 为 图 8 给 出 了 临界 温度 
平衡 态 。 所 以 , 极 化 初始 位 形 演化 为 有 序 状态 的 速 “附近 两 个 固定 温度 下 各 阶 矩 随 观测 时 间 的 变化 。 在 
度 比 随机 初始 位 形 快 。 显 然 , 初 始 位 形 偏离 平衡 状 。 “六 0 时 ,图 ga~d) 中 的 Cs 全 部 为 零 , 这 是 因为 初 
态 越 远 , 弛 耶 时 间 越 长 。 因 此 ,在 共 左 侧 , 随 机 位 始 位 形 是 随机 位 形 的 结果 。 当 初始 位 形 是 极 化 位 形 
的 平均 弛 豫 时 间 更 长 ,而 在 工 右 侧 , 极 化 位 形 的 平均 。 时 ,Ci 为 1,C,,4 为 零 ,如 图 8Ce~h) 所 示 。 


SS FP 


PRT MA EK, WTC) Piz. AT ARP 图 8(a) 和 (e) 中 ,C, 随 时 间 单 调 变 化 ,最 后 以 指 
化 的 最 大 影响 ,下 面 将 在 T<T 时 ,选用 随机 初始 位 。 数 接 近乎 衡 值 ,这 与 朗 之 万 方程 的 动力 学 演化 行为 
形 ;而 在 TT 时 ,选用 极 化 初始 位 形 。 一 致 。 表 明 Metropolis 算法 很 好 地 模拟 了 Ising 模型 


二 的 非 平衡 演化 。 图 8Ca~d)TT=0.99, 图 8Cb) 中 ,C, 先 
增加 后 减少 ,在 演化 过 程 中 形成 峰值 。 图 8Ce) 和 Cd) 

非 平衡 演化 对 观测 量 的 影响 可 以 从 观测 时 间 的 。 ”中 ,C, 和 C, 在 接近 稳定 信之 前 都 经 历 振 落 , 振 东 导 
角度 来 看 。 我 们 知道 弛 瑟 时 间 的 分 布 随 温 度 变 化 ， ”至 演化 过 程 中 的 符号 变化 。C, 和 C, 的 符号 可 以 是 
对 于 国定 温度 和 尺寸 的 系统 ,两 次 模拟 的 弛 殉 时 间 。 正 或 负 , 正 或 负 取 决 于 演化 时 间 。 


| Polarized Config. 


«Osh [ee T/T.= 1.01 


C,/ 10° 


1 2 
t/ 10° 


2 0 1 2 1 
t/10° t/ 10° t/10° 


8 在 7/7=0.99 (a~d) All 1.01 (e~b) 时 ,给 定 系 统 大 小 二 60 时 | 四 各 阶 矩 的 时 间 演 化 
Fig.8 The time evolution of cumulants of |m| at a given system size L=60 at 7/7,=0.99 (a~d) and 1.01 (e~h) 


Fl 8 Ce~h) T/T=1.01 , BT FE REIN TE AR HR C, tas 5R 8 Ce) RA FF SL, (LE 
类 似 于 图 8(a~d)。C, 在 接近 稳定 值 之 前 经 历 非 单 。 图 8(g) 与 (c) 相 反 , 符 号 首先 为 负 , 然 后 变 为 正 。 图 
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8() 中 的 CG 也 以 与 图 8(d) 类 似 的 符号 变化 , 先 为 
负 , 然 后 为 正 。C, 和 C, 的 符号 也 可 以 是 正 或 负 , 正 
或 负 取 决 于 观测 时 间 。 这 里 ,C, 和 C, 的 符号 变化 与 
动力 学 方程 的 解 一 致 2 1。 

但 是 ,要 注意 到 定量 上 工 两 侧 各 阶 矩 的 差异 非 
常 大 。 首 先 ,在 TIT=0.99 时 ,接近 C, 平 衡 预期 所 需 


混合 样本 提供 了 一 个 很 好 的 背景 信号 。 同 时 ,动力 
学 高 阶 算 能 有 效 观 测 RHIC BES vars 界 相关 的 
涨 落 。 


2) 对 于 有 限 系 统 尺 度 对 临界 观测 量 的 影响 , 首 
先 利 用 三 维 三 态 Potts 模型 计算 了 在 不 同 相 边界 区 
域 ,不 同系 统 尺 度 时 ,磁化 率 高 阶 矩 随 温 度 的 涨 落 ， 


的 时 间 约 为 4000, 而 在 T/T=1.01 时 ,所 需 时 间 小 于 
2 000。 差 异 超 过 两 倍 。 也 就 是 说 ,对 于 7< 工 的 系 
统 来 说 ,在 相同 的 时 间 内 达到 平衡 需要 更 长 的 时 间 。 
另外 ,C 和 C 在 7< 和 时 的 振荡 幅度 远大 于 7> 开 , 相 
差 大 约 是 两 个 数量 级 。 这 意味 着 在 高 于 临界 温度 
端 , 非 平衡 效应 的 影响 非常 小 ,可 以 忽略 不 计 ; 而 在 


低 于 临界 温度 端 , 非 平 衡 效 应 的 影响 非常 大 ,不 可 
4 结语 


我 们 系统 研究 了 相对 论 重 离子 碰撞 中 ,由 观测 
量 一 一 守恒 荷 高 阶 矩 确定 临界 涨 落 时 , 非 临 界 涨 落 、 
有 限 系统 尺度 和 有 限 演化 时 间 对 观测 结果 的 影响 。 


发 现 磁 化 率 二 阶 矩 在 一 级 相 变 、 二 级 相 变 临 界 点 以 
及 平 消 过 渡 区 域 都 具有 相似 的 峰值 结构 ,而且 ,峰值 
的 位 置 随 系 统 的 扩 度 变化 ;同时 ,磁化 率 四 阶 矩 在 这 
三 个 区 域 也 都 出 现 符号 改变 和 振荡 的 行为 。 表 明 在 
相 变 温度 附近 ,有 限 系统 高 阶 感应 率 中 会 出 现 非 单 
调 , 或 者 符号 改变 的 行为 ,这 些 行为 不 仅 发 生 在 二 级 
相 变 的 临界 温度 附近 ,也 会 发 生 在 一 级 相 变 和 平滑 
过 渡 区 域 。 而 且 , 对 于 有 限 系统 尺度 系统 ,其 临界 参 
数 也 不 能 简单 地 由 观测 量 峰值 的 位 置 决定 ,必须 由 
有 限 尺 度 的 标 度 规律 确定 。 
因此 ,提出 利用 有 限 尺度 标 度 律 所 表现 的 固定 
点 行为 能 更 精确 地 给 出 相 变 参数 ,并 给 出 了 对 固定 
点 行为 的 定量 描述 , 即 : 定 义 了 不 同 尺度 标 度 观测 量 
在 相 变 参数 附近 的 分 布 宽度 ,宽度 的 极 小 值 对 应 于 


1) 对 于 非 临界 涨 落 ,讨论 了 和 守恒 荷 高 阶 矩 测量 


司 定点 。 反 之 ,由 这 样 定义 的 宽度 在 相 变 参数 (温度 


所 需要 的 基本 统计 量 。 发 现 对 于 四 阶 和 矩 ,如果 按 每 
个 Nch 的 CBWC 方法 计算 ,总 统计 量 必须 大 于 
25 M。 建 议 用 0.1% 中 心 度 作 为 Nch bin 的 宽度 ,不 
仅 能 与 每 个 Nch 一 样 精确 代表 初始 体积 的 涨 落 ,而 
且 , 所 需 的 总 统计 量 在 1 M 量 级 就 可 以 保证 四 阶 矩 
测量 的 可 靠 性 。 对 于 六 阶 矩 ,总 统计 量 比 四 阶 矩 增 
加 一 个 量 级 (300 M) 都 不 够 ,0.1% 中 心 度 的 建议 对 
六 阶 矩 是 否 可 行 尚 需 进一步 研究 。 

对 于 末 态 粒子 数 有 限 所 导致 的 统计 涨 落 ,建议 
从 所 有 粒子 和 反 粒 子 独立 发 射 并 服从 泊 松 分 布 假设 
出 发 ,推导 了 三 种 守恒 荷 的 分 布 和 高 阶 窍 ,结果 都 和 
玻 尔 效 曼 近似 下 的 强 子 共振 气体 模型 所 得 到 的 结果 
一 致 。 通 过 将 泊 松 分 布 估计 的 统计 涨 落 与 RHIC 上 
产生 的 净 质 子 数 高 阶 矩 比较 发 现 , 实 验 测量 的 净 质 
子 数 高 阶 矩 对 碰撞 中 心 度 和 能 量 依赖 性 主要 是 统计 
涨 落 的 贡献 ,并 强调 了 统计 涨 落 不 可 忽略 ! 

最 后 ,给 出 了 用 池子 混合 事件 样本 的 方法 估计 
整体 和 系统 非 临界 涨 落 对 守恒 荷 高 阶 矩 的 影响 , 定 
义 动 力学 高 阶 矩 为 实验 数据 样本 所 测量 的 高 阶 矩 减 
去 对 应 混合 样本 的 高 阶 窍 。 通 过 利用 AMPT default 
模型 模拟 产生 19.6 GeV 金 - 金 磁 撞 的 样本 ,展示 了 池 
子 混合 样本 方法 的 有 效 性 。 结 果 表 明 : 所 得 到 
动力 学 高 阶 矩 不 仅 能 有 效 扣 除 统计 涨 落 的 影响 ,而 
ee bin 宽度 修正 和 探测 器 效率 修正 公 
式 一 样 有 效 扣 除 这 两 个 因素 对 高 阶 矩 测量 的 影响 。 
因此 ,对 于 由 全 局 和 系统 效应 所 导致 的 非 临 界 涨 落 ， 


和 相 变 指数 率 ) 平 面 的 等 高 图 ,可 以 精确 决定 相 变 
参数 。 

以 三 维 三 态 Potts 模型 的 平均 磁化 强度 为 例 ,名 
出 了 平均 磁化 强度 在 相 边 界 二 个 不 同 区 域 附近 , 标 
度 曲 线 在 相 变 参数 平面 上 的 等 高 图 ,由 图 的 分 布 极 
值 位 置 ,以 及 分 布 特征 ,不 仅 可 以 精确 确定 相 变 参 
数 , 而 且 可 以 判别 相 变 级 数 。 

当然 ,将 固定 点 方法 应 用 到 相对 论 重 离子 碰撞 
中 还 有 很 多 待 解决 的 问题 ,首先 是 必须 选择 跟 磁 化 
率 类 似 的 观测 量 , 守 恒 荷 高 阶 和 矩 可 能 是 一 个 非常 好 
的 选择 。 然 后 是 如 何 将 实验 中 的 碰撞 能 量 、 中 心 度 
等 与 系统 的 温度 、 尺 度 对 应 。 

3) 对 于 有 限 演化 时 间 导 致 某 些 系统 可 能 没有 达 
到 热力 学 的 平衡 态 , 即 非 平衡 演化 过 程 的 影响 ,提出 
利用 Metropolis 算法 模拟 三 维 伊 辛 模型 的 非 阅 衡 演 
化 ,并 定义 了 系统 的 平均 弛 豫 时 间 ,通过 模拟 零 外 磁 
场 下 三 维 伊 辛 模 型 的 非 平衡 演化 ,结果 显示 : 伊 辛 模 
型 中 序 参量 的 非 平衡 演化 确实 以 指数 方式 接近 其 稳 
定 值 ,与 朗 之 万 方程 的 解 相同 。 临 界 温度 下 的 平均 
弛 豫 时 间 随 系统 尺度 的 z 次 过 发 散 , 与 动力 学 方程 
中 的 弛 豫 时 间 相 同 ,表明 所 用 模拟 方法 与 解 动力 学 
方程 一 样 有 效 , 所 定义 的 平均 弛 鸳 时 间 正 确 描 述 了 
系统 动力 学 时 间 演 化 的 特征 。 

同时 ,通过 模拟 方法 研究 了 非 平衡 演化 的 特征 ， 
发 现 平均 弛 驳 时 间 取 决 于 系统 尺度 和 温度 。7>ZT 时 
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的 平均 弛 更 时 间 非 常 短 , 几乎 与 初始 位 形 无 天; 在 7 
<T 时 ,当初 始 位 形 远 离 平 衡 状态 时 ,平均 弛 豫 时 间 
比 7>7T 长 很 多 ,在 一 级 相 变 线 上 ,系统 更 难 达 到 平 
衡 。 因 此 ,在 一 级 相 变 区 域 , 非 平 衡 因 素 的 影响 不 可 
忽略 。 但 是 ,在 平滑 过 渡 区 域 , 非 平 衡 因 素 的 影响 可 
以 忽略 不 计 。 

通过 研究 非 平 衡 演 化 对 观测 量 的 影响 ,给 出 了 
序 参 量 各 阶 矩 的 时 间 演 化 结果 。 结 果 发 现 非 平衡 过 
程 对 高 阶 矩 的 影响 比 低 阶 矩 大 ,并 且 ,C, 和 C 都 显 
示 振 荡 ,其 值 可 能 取 正 或 负 , 取 决 于 观察 时 间 ,这 与 
动力 学 理论 计算 的 结果 一 致 。 

这 些 非 平衡 演化 的 定性 特征 趋势 不 依赖 于 系 
统 , 可 以 直接 作为 相对 论 重 离子 磁 撞 实验 测量 的 参 
考 , 非 平衡 演化 对 确定 QCD 相 边 界 的 影响 在 平滑 过 
渡 区 域 几乎 可 以 忽略 不 计 , 但 是 ,在 一 阶 相 变 线 上 ， 
非 平 衡 演化 的 影响 不 可 忽略。 另外 ,由 于 非 平衡 态 
的 影响 , 序 参 量 的 三 阶 矩 的 符号 可 能 出 现 负 值 , 与 
RHIC/STAR 的 测量 结果 一 致 。 

本 文 综述 了 近 15 年 来 ,我 们 在 通过 相对 论 重 离 
子 人 碰撞 确定 QCD 相 边 界 方面 的 主要 工作 。 由 此 可 
见 , 最 终 确定 QCD 相 边 界 仍然 任重道远 ,还 有 很 多 
问题 需要 进一步 的 研究 。 关 于 间 欣 理论 和 实验 方面 
的 工作 可 以 参考 综述 文章 中 。 
作者 贡献 声明 吴 元 芳 主 要 负责 本 文 的 撰写 ,总 体 
设计 和 指导 了 大 部 分 的 具体 工作 ; 李 笑 冰 8$3 工作 的 
主要 完成 人 ,主要 负责 该 论文 的 编辑 .审核 工作 ; 陈 
丽 珠 $1.1 工 作 的 主要 完成 人 ,完成 81.2 的 部 分 工作 ， 
与 $1.3 以 及 8$2 部 分 工作 的 讨论 ,并 完成 了 部 分 黄 
基 工 作 ; 李 治 明 参与 部 分 工作 的 讨论 ,提出 建议 ; 许 
明 梅 $3 部 分 工作 的 主要 撰 稿 、 指 导 , 参 与 其 他 部 分 
工作 的 讨论 ,提出 建议 ; 潘 雪 $1.2、82.1 工作 的 主要 
完成 人 ,并 完成 了 8§2.2 的 奠基 工作 ; 张 几 8§1.3 工 作 的 
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